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ABSTRACT 



Context, Spatially resolved observations of circumstellar disks at millimeter wavelengths allow detailed comparison with theoretical 
models for the radial and vertical distribution of the material. 

Aims. We investigate the physical structure of the gas component of the disk around the pre-main-sequence star HD169142 and test 
the disk model derived from the spectral energy distribution. 

Methods. The I3 CO and C 18 7=2-1 line emission is observed from the disk with 1"4 resolution using the Submillimeter Array. We 
adopt the disk physical structure derived from a model which fits the spectral energy distribution of HD169142. We obtain the full 
three-dimensional information on the CO emission with the aid of a molecular excitation and radiative transfer code. This information 
is used for the analysis of our observations and previous 12 CO J =2-1 and 1.3 mm continuum data. 

Results. The spatially resolved 13 CO and C ls O emission shows a Keplerian velocity pattern The disk is seen at an inclination close 
to 13° from face-on. We conclude that the regions traced by different CO isotopologues are distinct in terms of their vertical location 
within the disk, their temperature and their column densities. With the given disk structure, we find that freeze-out is not efficient 
enough to remove a significant amount of CO from gas phase. Both observed lines match the model prediction both in flux and in the 
spatial structure of the emission. Therefore we use our data to derive the 13 CO and C ls O mass and consequently the 12 CO mass using 
standard isotopic ratios. 

Conclusions. We constrain the total disk gas mass to (0.6-3. 0)xlO~ 2 M Q . Adopting a maximum dust opacity of 2 cm 2 g~ x we derive 
a minimum dust mass of 2.16xl0~ 4 M Q from the fit to the 1.3 mm data. Comparison of the derived gas and dust mass shows that 
the gas to dust mass ratio of 100 is only possible under the assumption of a dust opacity of 2 cm 2 /g~' and 12 CO abundance of 10~ 4 
with respect to H 2 . However, our data are also compatible with a gas to dust ratio of 25, with a dust opacity of 1 cnr/g -1 and l2 CO 
abundance of 2xl0 -4 . 
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1. Introduction 

Although the presence of molecular gas in disks around interme- 
diate mass pre-main-sequence (Herbig Ae) stars was reported a 
decade ago (Mannings & Sargent) 1997I) . the research in this field 
has focused more on their less massive counterparts (T Tauri 
stars). The disks ar ound T Tauri stars hav e masses ranging from 
0.001 to 0.1 M Q dBeckwith et al.l[T990h . usually derived from 
millimeter continuum fluxes assuming a gas to dust mass ratio 
of 100, as in molecular clouds. Their outer radii are constrained 
by m olecular line observations and are typically a few hundred 
AU (ISimon et al.ll2000t iThi et al.ll2001l and references therein). 
Due to the low luminosity of the central star (0.5 to 1 L ), these 
disks are relatively cold (less than 20 K beyond 100 AU from the 
star) causing a significant depletion of the CO in the midplane of 
the outer disk. On the other hand, the Herbig Ae stars are about 
ten times more luminous than T Tauri stars, and consequently 
their disks are warmer. This allows the easiest to detect and the 
most commonly used gas tracer, carbon monoxide, to remain in 
the gas phase even in the disk midplane. Observations of CO and 
its isotopologues toward Herbig Ae stars are therefore expected 
to be more powerful probes of the full disk structure. Only a 
few Herbig Ae disks have been studied thoroughly via spatially 
resolved observations of molecular line emis sion that include s 
the opticall y thin CO isotopo logues: AB Aur (Pietu et al. 2005), 
MWC480 (Pi etu etal.ll20 07). HD 163296 (Isellaet al. 2007b- 



2. HD169142 

The object of our study, HD169142, is a 2.0 M Herbig Ae 
star of spectral type A5Ve sur rounded by a gas-rich circum- 



stellar disk loca ted at 145 pc (ISvlvester et al . 1996). With an 



age of 6^3 Myr dGradv et al.ll2007l) and its spectral energy dis 
tribution marked by infra red excess and the lack of silicate 
features dDent et all 120061) . HD169142 is an example of an 
advanced pre-main-sequence evolutionary stage. Unlike most 
of the Herbig Ae/Be stars, it shows no e vidence of proxim - 
ity to a cloud or extended molecular gas (Me eus et al.l 12.00 lb . 
Observations of molecular gas in this disk are therefore eas- 
ier to interpret. Howev er, HD 169 142 is no t completely isolated 
from other young stars: Grad y et all d2007l) find three coeval pre- 
main-sequence stars within a projected separation of 1160 AU. 
The closest companion is located at 9'.'3 separation and may 
form a binary system with HD169142. Near-infrared polarisa- 
tion imag es show that th e dust in the disk extends to at least 
217 AU (Kuhn et al. 200 J). More rec ent submillimeter observa- 
tions (Dent et al. 2005; Raman et al. 2006) show bright and nar- 
row CO lines. iRaman et al.l (|2006) spatially resolve the disk and 
find a fit to the CO 7=2-1 line and 1.3 mm continuum observa- 
tions by adopting a flared accretion disk model with a 235 AU ra- 
dius and a 13° inclination from face on. Observations at optical, 
IR, and (sub)millimeter wavelengths allow ed modelling of the 
disk's spectral energy distribution (SED) (MalfaitetaL|[l998; 
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Fig. 1. Integrated intensity (contours) and first moment maps (colour scale) of 12 CO J=2-l (left panel, from Ram an et al.l (H006)), 
13 CO J=2-l (middle panel) and C 18 7=2-1 line (right panel). Contours are 1, 2, 3,...x200 mJy beam 1 km s" 1 for 12 CO and 13 CO, 
and 1, 2, 3,...xl00 mJy beam -1 km s~' for C 18 0. The integrated intensity and first moment maps are obtained over a velocity range 
of 5.6-8.4 km s 1 . The data were clipped at 0.7, 0.5, and 0.35 Jy beam 1 for 12 CO, 13 CO and C 18 Q, respectively. 



iDominik et al1l2003t iDent et al.1120061) . iMalfait et all (1 19981) fit- 
ted the near-infrared and far-infrared excess of HD169142 by 
two disk components: an inner disk extending from 0.5 AU to 
1 AU with a density exponent of 2.0 an d the outer disk from 
28 AU with a flatter density distribution. IDominik et al.l 12003) 
adopt a low inclination of 8°, outer radius of 100 AU and sur- 
face density exponent p-2 to fit the SED, and therefore derive 
a d isk mass of 0.1 M H . A more detailed SED modeling is done 
by IDent et all d2006l) where both the SED and resolved 7 mm 
continuum emission w ere fitted using an accretion disk model 
dD'Alessio et al.1120051) corresponding to a 10 Myr old A2 spec- 
tral type star. They adopt an inclination of 30° and outer radius 
of 300 AU, and derive a disk mass of 4xl0~ 2 Mn. lGrady et alJ 
d2007l) fit the SED and NICMOS image at 1.1 ^m with a model 
consisting of two distinct disk components - the inner disk from 
0.15 to 5 AU radius and the outer disk extending from 44 to 
230 AU. It is important to stress that all above mass estimates 
of the disk around HD169142 are based solely on the observed 
dust emission, and not gas. 

This paper presents resolved interferometric observations of 
the 13 CO and C 18 7=2-1 lines from HD169142. The observa- 
tions and results are shown in Sections 3 and 4. Section 4 intro- 
duces the disk model we a dopt dD'Alessio et al.l2005l : lDent et alJ 
2006: lRaman et al.| [2006) and our fit to the 1.3 mm data provid- 
ing a measure of the minimum disk dust mass. We fit the 13 CO 
and C 18 emission deriving the corresponding 13 CO and C ls O 
mass in the disk, and place constraints on the total gas mass of 
the disk. We discuss the implications of our findings for the gas 
to dust ratio in the disk and the possible presence of substructure 
in the disk midplane. Section 5 summarizes our findings. 

3. Observations and results 

The observations of HD169142 were carried out with the 
Submillimeter Arra}{3 (SMA) on 2005 April 19, simultane- 
ous with the obser vations of 12 CO J=2-l line presented in 
Ram an et alJ |2006). A more detailed description of the obser- 
vations and of the calibration procedure is given there. The cor- 

1 The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 



relator provided 2 GHz of bandwidth in each sideband and was 
configured to include the 13 CO J=2-l line at 220.3986765 GHz 
and the C ls O J=2-l fine at 219.5603568 GHz in the lower side- 
band in a 104 MHz wide spectral band with channel spacing of 
0.2 MHz (-0.26 km s" 1 ). 

The data reduction and i mage analysis we re done with the 
Miriad data reduction tools dSault et alj|1995l) . The (u,v) data 
were Fourier transformed using natural weighting. The resulting 
synthesized beam size is 1'.'4 x l'.'O (PA=26°). The rms of the 
line images is 180 mJy beam -1 per channel or 4.6 K (4.8 K for 
12 CO). 

Emission of 13 CO and C 18 7=2-1 was detected from the 
HD169142 circumstellar disk. Figure Q] shows the intensity 
weighted velocity maps with overlaid integrated intensity con- 
tours for bo th lines as well as the previously published 12 CO 
7=2-1 line dRaman et al.ll2006l) . The intensity integrated over 
the velocity range from 5.6 to 8.4 km s , in which the line 
emission is fully contained, and over the central 4"x4" region 
is 12.1, 6.5 and 2.7 Jy km s 1 for 12 CO, 13 CO and C 18 7=2- 
1 line, respectively. All three lines follow a similar velocity 
pattern, interpreted as a clear indication of Keplerian rotation 
around a 2 M star of a disk seen at 13° inclination (Raman 
et al., 2006; this work). Figure [2] shows the 12 CO, 13 CO and 
C 18 7=2-1 line spectra summed over a 4"x4" region centered 
at the HD169142 position. The profiles are relatively symmet- 
ric, double peaked and centered at 7. 1 +0.2 km s _1 and reflecting 
the underlying rot ation pattern. The 12 CO and 1.3 mm results 
were presented by Ra man et al.l ([2006). They report a 1.3 mm 
line flux of 169+5 mJy. In Fig.[3]the emision stucture is shown 
over a range of velocities corresponding to the 13 CO 7=2-1 line 
(top panel). The emission extends to 2" (270 AU) from the star 
at a 2<r level and shows a Keplerian velocity pattern. Figure [4] 
presents the channel maps of the C ls O 7=2-1 line (top panel) 
with a kinematic structure similar to that of the 13 CO 7=2-1 line 
seen in Fig. [3] The 2cr level emission reaches 1'.'5 (220 AU) from 
the star. At 7.4 km s _1 a localised emission peak of 0.94+0.35 Jy 
is seen l'.'O north from the star. 



4. Discussion 

iRaman et all d2006) show that the model that fits the SED of 
HD169142, can also be used to fit the structure and intensity of 



O. Panic et al.: Gas and dust mass in the disk around the Herbig Ae star HD169142 



3 



~3 



X 

3 



I I I I I I I I I I I I I I I I I I I I I I I 



I 



,2 CO (2-1) 



,3 C0 (2-1) 



1 1 

10 
9 
8 
7 
6 
5 

u 4 
3 
2 
1 


-1 

2 4 6 8 10 12 
V LSR (km s- 1 ) 

Fig. 2. Spectra of C lg O J=2-l ( bottom), 13 CO J=2- l (middle) 
and 12 CO 7=2-1 line (top, from Ra man et alj d2006l) ) summed 
over the central 4"x4" toward HD169142. The 13 CO and 12 CO 
fluxes are shifted vertically by 2 and 5 Jy respectively. Black 
lines show the observed spectra and the red lines sh ow the model 
fit fou nd in Section 4.3. 1 (for 12 CO the model from Raman et al. 
d2006l) was used). 



the resolved 12 CO emission. They note that some weak residual 
emission is still present after subtracting the model from the ob- 
servational data, that may correspond to real substructure within 
the disk. We analyse 13 CO and C 18 data, which are expected to 
probe different depths in the disk due to their lower opacities. 

4.1. Adopted disk model 

iDent et all d2006l) fit the spectral energy distribution (SED) 
of HD169142 using an accretion disk model from the 
ID'Ale ssio et al. (2005) database with an accretion rate of 
10 -8 M yr~', an outer radius of 300 AU, and a 30° inclination. 
The age of the central star in this model is 10 Myr. They adopt an 




A2 spectral type for the central star because it provides a slightly 
better fit than the A6 spectral type model does to the SED in 
the range of 2-200 fim dominated by the warm dust emission. 
However, there is no significant difference between the SEDs 
corresponding to these two models. Beyond 44 AU from the 
star, the model is qualitatively similar to the mod el of the outer 
disk of HD169142 described in lGradv et alld2007l) . The detailed 
model structure, including plots of surface density distribution, 
temperature at different scale heights and optical thickness of the 
disk, is provided in the online database of accretion disk models 
(ID'Alessio et alJl2005lFI . 

In order to match thei r resolved submillimeter ob servations 
of CO , lRaman et all d2006l) modify the model used in lDent et alj 
(2006) to have a disk radius of 235 AU. Assuming the mass of 
the central star to be 2 M Q , they derive an inclination of 13+1° 
from face-on and use it in the model. These modifications are 
not expected to affect the quality of the SED fit since the SED 
of HD169142 alone does not constrain well the outer radius and 
inclination of the disk. They show that this SED model can be 
used to fit the structure and intensity of the resolved 12 CO emis- 
sion. However, the SED does not provide reliable constraints of 
the disk dust mass, which in the model is set to 2xl0~ 4 M . This 
parameter depends only on the optically thin thermal emission 
from the disk in the millimeter wavelength region. 

The density in the model is a parameter that can be scaled 
by a factor of few, without significantly affecting the computed 
disk temperature structure. The resulting mass only affects the 
millimeter flux in the modeled SED, i.e., the 1.3 mm data point- 
In sec tion 4.2 we re-analyse the 1.3 mm data from Rama net alj 
(2006) using a radiative transfer code to calculate the emission 
from the model and better constrain the dust mass of the disk. 

In the analysis of our 13 CO and C ls O line emission (sec- 
tion 4.3) we use this improved disk model, which already en- 
compasses a range of observational constraints. We assume that 
the gas follows the dust distribution prescribed by the model and 
that their temperatures are the same. The micro turbulent FWHM 
of the lines is set to 0.16 km s throughout the disk. One con- 
clusion that we can already draw from the disk model is the frac- 
tion of disk mass where freeze-out of CO is likely to be efficient, 
which is believed to happen at temperatures below 20 K. The 
mass weighted temperature in the model is 35 K, which indicates 
that the bulk of the disk material is at temperatures above 20 K. 
The temperature falls below 20 K only in the midplane region of 
the disk beyond 150 AU from the star and there the abundance of 
the gas phase CO and its isotopologues is expected to be heavily 
decreased due to freeze-out of these molecules onto dust grains. 
This region contains only 8% of disk mass and therefore any 
mass estimate based on 13 CO and C 18 emission may represent 
92-100% of the true disk mass. 

In all previously published work on HD169142, the disk gas 
mass was inferred indirectly from the derived disk dust mass, 
adopting a standard gas to dust ratio. In our analysis of the 
13 CO and C 18 emission we argue that it is reasonable to as- 
sume a standard 12 CO abundance in the case of HD169142, be- 
cause freeze-out affects less than 8% of the mass in this disk. 
Therefore we investigate the constraints that can be obtained 
from our molecular line data on the disk gas mass. 

4.2. Dust continuum emission 

From the observed 1.3 mm continuum flux of 169 ±5 mJy, 
Ram an et alj ([2006) derive a dust mass of 2xl0~ 4 M Q assum- 
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Fig. 3. - Top panel: Channel maps of the observed 13 CO 7=2-1 emission at a range of velocities over which the line is detected. 
The lower left corner of this panel shows the size and position angle of the synthesized beam. - Middle panel: 13 CO 7=2-1 chanel 
maps from the best-fit disk model from Section 4.3. 1 . - Lower panel: Channel maps showing the residual emission after subtraction 
of the best-fit model from the data. The contour levels are -1,1,2,3,4x360 mJy beam -1 (2 sigma) in all panels. 
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Fig. 4. - Top panel: Channel maps of the observed C ls O 7=2-1 emission at a range of velocities over which the line is detected. 
The lower left corner of this panel shows the size and position angle of the synthesized beam. - Middle panel: C ls O 7=2-1 channel 
maps from the best-fit disk model from Section 4.3. 1 . - Lower panel: Channel maps showing the residual emission after subtraction 
of the best-fit model from the data. The contour levels are -1,1,2,3,4,5x360 mJy beam -1 (2 sigma) in all panels. 



ing an emissivity k^, us '-2 cm 2 g du 1 st and a single temperature of 
30 K. We reanalyse the 1.3 mm emission using the disk tem- 
perature and density structure from the model described in sec- 
tion 4.1, where we vary the dust density (i.e., the dust mass) 
to fit the observations. A large uncertainty on the mass deter- 
mination from the 1.3 mm continuum flux is due to the dust 
opacity, which is not well d etermined in circumstellar disks. 
Ossenkop f & Henningl (l994) study the effects of dust coagu- 
lation, and ice coverage of grains on their opacity in protostel- 
lar cores and suggest 1.3 mm opacity of 1 cm 2 g du ' st for very 
dense (n>10 7 cm 4 ) regions, and up to 5 times more if the 
grains are not covered in ice mantles. The latter is thus only 
valid fo r regions ab o ve wa ter ice desorption temperature (80- 
100 K, iFraser et ail (120011)1 . In their study of the dust opac- 
ity in circumstellar disks Draind (120061) find that astrosilicate 
and pirolysed cellulose at 600 °C are materials representative 



of the dust properties which may be expected in circumstellar 
disks. They explore a large range of grain size distributions and 
find that these materials have mm-wavelength opacities which 
are close to the observational constraints on dust opacity from 
ex tinction studies of the diffuse interstellar medium presented 
in IWeingartner & Draine] d2001l) . The resulting opacity could 
be anywhere between 0. 1 and 2 cm 2 g du 1 st at the wavelength of 
1.3 mm, depending on the adopted grain size. We use this in- 
formation to calculate the minimum amount of dust needed to 
produce the observed flux by adopting the opacity (and emissiv- 
ity) of 2 cm 2 g du 1 st . The SED modelling of the disk emission at 
shorter wavelengths is not affected by our assumptions of dust 
opacity at 1.3 mm. The dust thermal emission is dominated by 
the cold disk midplane at large radii, likely to contain settled 
and/or grown dust particles whose optical properties may differ 
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from those of the small dust at disk surface, responsible for the 
optically thick near-infrared and infrared emission of the disk. 

We do the full modeling of th e interferometric visibilities 
using the radiative transfer code of Hogerheiid e & van der Takl 
(2000) (RATRAN). We perform a x 2 minimisation to fit the 
1 .3 mm visibilities by varying the disk dust mass. The best match 
is shown in Fig. [5] top-left panel, resulting in an estimate of the 
dust mass of the disk, given by 



Mdust 



2 cm 2 
= 2.16x 10~ 2 ^ 

^1.3 mm 



2 CITIES 

=xM . 



(1) 



Assuming that the adopted dust opacity of 2 cm 2 g du 1 st is the 
maximum value, our estimate presents the lower limit on the 
dust mass in HD169142: M dust >2.16x 1O~ 4 M . A slightly 
higher maximum dust opacity of 7 cm 2 g du I st is obtained for py- 
rolised cellulose at 800 °C but this value exceeds greatl y the 
observational contraints from lWeingartner & Dra ine (2001 ]) an d 
although not impossible, is unlikely in disks (Draine 2006). 
However we cannot exclude the possibility that other dust prop- 
erties like cristallinity, porosity and shape, which are not well 
understood, may enhance dust submillimeter opacity beyond 
2 cm 2 g, 1 by a factor of a few. 



4.3. Molecular line emission 

Due to their low abundance relative to 12 CO, the 13 CO and C ls O 
molecules emit lines that are comparatively less optically thick. 
As a result of this, they saturate at different heights in the disk 
allowing us to use them as a probe of the physical conditions 
in different disk lay ers and investigate the disk vertical structure 
dDartois et al.l l2003). Furthermore, the C I8 line emission is of- 
ten optically thin in the outer disk region, where the bulk of disk 
mass is located, which makes it a valuable probe of the disk gas 
mass. 

The commonly employed route for the analysis of the CO 
isotopologue and dust continuum emission from circumstellar 
disks is to convert the continuum flux to disk mass adopting a 
certain dust opacity, usually 0.02 cm 2 g ga ' s (iBeckwith et al.ll 990). 
under the assumption of the canonical gas to dust mass ratio of 
100. In this way, the molecular line emission is fit for abundance, 
where any discrepancy with respect to the canonical abundances 
and isotopic ratios of observed molecules is explained by a de- 
pletion factor, indicating the fraction of the CO gas affected by 
the freeze-out onto dust grains or selective photodissociation. 
This approach is used in lRamanet al. (2006) where the analysis 
of the 12 CO J=2-l line emission and 1.3 mm continuum from 
HD169142 is done for a disk mass of 2xl0~ 2 M Q and 12 CO was 
found to be depleted by a factor four, having an abundance of 
2.5xl0~ 5 . This corresponds to 6xl0~ 6 M Q of 12 CO molecules 
present in gas phase in the disk. In the present work, we argue 
that the 12 CO abundance is not significantly affected by freeze- 
out or selective photodissociation and we chose to follow a dif- 
ferent route, in which we derive a rough gas mass estimate from 
the molecular line observations taking into account all related 
uncertainties. 

In order to interpret our 13 CO and C 18 observations, we 
adopt the disk model described in detail in section 4. 1 . We set 
the dust mass in the model to 2.16xl0~ 4 M ; the value derived 
as the minimum dust mass for HD169142 in section 4.2. We as- 
sume that the molecular gas follows the dust temperature and 
density structure and we vary the 13 CO and C 18 gas masses to 
best fit the observations. We calculate the 13 CO and C 18 7=2- 
1 line emission from the model using the RATRAN molecular 



excitation and radiative transfer code. The model visibilities are 
generated with the miriad task uvmodel, and the comparison to 
the observed visibilities is done directly. We use assumptions 
about the molecular abundances of the 12 CO and isotopic ratios 
to make a conversion to the total gas mass (1.2xMh,)- We ex- 
plore a range of values of the gas to dust mass ratio fg/j. 

4.3.1. Disk structure probed by 13 CO and C 18 

We find fits to the observed visibilities of 13 CO and C 18 line for 
Mi3 CO =2.9xlO~ 7 M and M c is o =4.6xl0~ 8 M , shown in Fig. [5] 
bottom panels. For comparison, the visibilities of the 12 CO J=2- 
1 line are also shown in Fi g. [5] top-right panel, with the best 
fit obtained by Ram an et al. (2006) for Mi2 CO =6xlO- 6 M con- 
tained in the gas phase. In order to test our fit further, we invert 
and deconvolve the visibilities using the miriad reduction pack- 
age and compare the observed (top panels) and modeled (mid- 
dle panels) velocity channel maps in Figs. [3] and [4] The struc- 
ture and intensity of the emission is found to be well matched. 
Figures [3] and [4] (lower panels) show the residual emission re- 
sulting from subtraction of the modeled from the observed data. 
For 13 CO the residuals appear randomly distributed and do not 
exceed the noise level. Significant residuals are seen in C ls O at 
7.4 km s -1 about 0'.'5 north from the star. The C ls O data (Fig. |4| 
have a marginally significant localised excess l'.'O north of the 
star seen at 3 cr level in the residual emission. This excess emis- 
sion corresponds to the redshifted peak of the spectral line at 
Vlsr-7-4 km s _1 . It is marginally resolved at 1'.'5 x 2'.'5 and con- 
tains about 20% of C 18 line flux. Similar, but not spatially co- 
incident residual emission is also seen in 13 CO (Figure [3] lower 
panel) near the noise levels, but not in 12 CO. A cause of this fea- 
ture may lie in disk midplane asymmetry. Any connection with 
the reported companion at 9'.' 3 is unlikely, but a hypothesized 
body within the inner disk is an interesting prospect ( Grad vet al.l 
2007). Nevertheless, an instrumental or calibration artefact may 
have affected the 7.4 km s channel resulting in a slight increase 
of flux in that channel in the lower sideband, only noticeable in 
the weak C 18 line. 

In this section we will use our 13 CO and C 18 results to make 
an estimate of the total gas mass M gas in the disk, related to the 
12 CO gas mass Mi2 C0 through 



1 r, 

1.2 x — x 

mi2 C0 



1 



12 CO] 



xM 



2 CO- 



(2) 



where m H , and mi2 C0 are the masses of the H2 and 12 CO 
molecule, and [ 12 CO] is the abundance of 12 CO with respect 
to H2. The factor of 1.2 accounts for one fifth of the gas con- 
tained in helium. We can use the derived masses of the 13 CO 
and C 18 gas to calculate the mass of the 12 CO gas. We adopt 
the isotopic ratios of 12 C/ 13 C=77±7 and 16 O/ 18 O=560+25 from 
IWilson&Roodl(fT99ll) . 

M,2 C0 = ™^ x ^ x M c , 80 = (2.4 ± 1 .0) x 10- 5 M o (3) 

M.2 C0 = ™^ x x Mnco = (2.2 ± 0.6) x 1O- 5 M (4) 
mu C0 [ 13 C] 

The values obtained here are given with errors which are domi- 
nated by the rms of our data and also include the errors on iso- 
topic ratios. If we take into account that 8% of the 12 CO and its 
isotopologues in HD169142 may be depleted onto dust grains 
and therefore not contributing to the emission observed, then the 
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Fig. 5. Top-left: Vector-averaged continuum flux as a function of projected baseline length (black dots). The error bars represent 
the variance within each annular average. The zero-signal expectation value is shown by the dashed histogram. The red dotted line 
shows the best-fit model for a dust emissivity of 2.0 cm 2 glL with the disk mass of 2.16xl0~ 4 M Q . Top-right: Same for the 12 CO 
emission from Ra man et al.l d2006l) . integ rated over the width of «2.7 km s _1 (10 spectral channels) centered on each line. The red 
dotted line shows the model presented bv lRaman et alj ((2006). Bottom-left and bottom-right: Same for our 13 CO and C ls O J=2-l 
data, respectively. The red dotted lines show the best-fit model from Section 4.3.1. 



full i2 CO mass range is (1.4-3.7)xl0 -5 M Q . The fact that both 
datasets roughly agree on the 12 CO mass provides an additional 
argument for the validity of the model . How ever, the fit to the 
12 CO J=2— 1 data done in iRaman et al.l d2006l) corresponds to a 
mass of gas phase 12 CO of 6xl0~ 6 M Q (mentioned in section 
4.3), two to six times lower than the range derived here. The 
12 CO emission is optically thick and therefore it does not trace 
the entire disk but originates from a layer much higher i n the disk 
than either the 13 CO or C ls O lines dDartois et alj2003l) . The dis- 
crepancy between the 12 CO result and that of the isotopologues 
indicates that the adopted disk model may not provide the best 
description in these high layers. We will discuss this in more de- 
tail in section 4.5 where the effects of temperature and turbulence 
on the line profile are discussed. 

Our model calculations allow us to identify which regions 
of the disk contribute most to the observed lines. We follow the 
line of sight integration step by step from the observer into the 
disk and find the location where 90% of the integrated intensity 
is reached for a face on orientation of the disk. We show these 
locations in Fig.|6]in a cross-section of the disk showing the un- 
derlying density structure. We choose to use the intensity inte- 
grated over a range of wavelengths in which the line is emitted in 
order to better represent the region corresponding to the physical 
quantity that is observed, i.e., the line flux. The use of surfaces 



with equal opacity (often used t—3) may be misleading, because 
the opacity at the line centre is not necessarily consistent with the 
opacity in the line wings, and the amount of flux still produced 
in the line wings beyond the t=3 surface may contribute signif- 
icantly to the line flux in some species more than in the others. 
The calculation was done for 12 CO, 13 CO and C ls O J=2-l line. 
The discontinuous shape of the surface above which 90% of the 
integrated intensity in these lines is emitted is due to the discrete 
cells in the grid of our model that are used to sample the den- 
sity and temperature structure. As expected, the less optically 
thick species trace regions deeper in the disk with C ls O trac- 
ing the disk midplane while 12 CO traces the disk surface layer. 
We analyse the physical properties of the emitting regions in the 
outer disk, beyond 100 AU from the star (scales resolved by our 
SMA observations). For the two outermost radial positions in 
our model grid, we calculate the column of 12 CO gas contained 
in the region where 90% of the emission is coming from, which 
is located above the vertical position shown in Figj6j for each 
molecular line. These values are shown in Table Qflj. Also, the 
average temperature of these gas columns is calculated, taking 
into account the temperature weighted by the column density 



3 The calculations f or 12 CO J=2-l are done for the model presented 
in Ra man et "all (l2006h 
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Fig. 6. Density structure of the best-fit disk model shown as colour scale. The disk is viewed by us from almost straight above. The 
three dotted lines show the heights above which 90% of the integrated line emission of 12 CO, 13 CO and C 18 7=2-1 originates. For 
the C 18 emission this line lies below the disk midplane. 



(i.e., mass) for each cell located above the vertical positions from 
Fig-El The C ls O J=2-l line generally traces three times as much 
CO gas as traced by the 13 CO J=2-l line in the outer disk while 
the 12 CO J=2-l line traces only a small fraction of the gas, a few 
percent of the column traced by C ls O. It is clearly seen that the 
C ls O traces deep into the disk, to the height of several AU at the 
far side of the midplane, showing that it is very sensitive to the 
total gas mass. The 90% of the C 18 J=2-l line emission traces 
around 60% of the total gas column. 

In section 4.1 we already concluded from the disk model 
that the CO freeze-out could only be efficient in the cold outer 
disk midplane region which contains 8% of the total disk mass. 
This very small fraction is due to the luminosity of the Herbig 
Ae star HD169142 which is higher than the typical luminosity 
of T Tauri stars. In the latter, the circumstellar disks are colder 
and the CO freeze-out presents a major obstacle to determine 
gas mass based on obse rvations of CO o r its isotopologues. As 
a caveat, we note that iPietu et ail d2006l) find somewhat lower 
dust temperature in the disk around the Herbig Ae star MWC 
480 (10-25 K at 100 AU). Since we have no indications to the 
contrary, we adopt the dust temperature that SED disk model- 
ing suggests. If observed with a sufficiently high resolution and 
sensitivity, disks where CO freeze-out is efficient are expected to 
have smaller radii derived from CO isotopologues compared to 
the outer radius derived from the 12 CO and from the dust emis- 
sion. Although our observations are limited by the beam size 
(l'.'4*;200 AU) we have tested a range of different outer radii 



and found that, based on our molecular line data, the outer ra- 
dius should not be very different from 200 AU for either of the 
observed lines. We also do not find any significant discrepancy 
between the molecular line data and dust continuum that could 
indicate different outer radii for gas and dust. 

Another factor that can influence the 12 CO abundance is pho- 
todissociation. The effect of photodissociation via the stellar UV 
radiation on CO abu ndan ce in circumstellar disks is studied in 
Ivan Zadelhoffl d2002l) and Ivan Zadelhoffet all d2003l) . They find 
that the CO abundance decrease due to photodissociation is sig- 
nificant only in the upper disk lay ers and that the CO c olumn 
density is not much affected by it. Jonkheidet al. (2004) show 
that the stellar UV field is only able to heat 10% of the gas 
mass in the disk to considerably high temperatures, in the sur- 
face disk layer. The effect of photodissociation on the CO con- 
tent of HD169142 is expected to be small, due to the significant 
dust mass of the disk. Since we adopt a maximum emissivity and 
therefore minimal dust mass, the conclusion that photoioniza- 
tion can be neglected is robust. In addition, the C 18 emission 
that we use to derive the gas mass, arises predominantly from 
the outer disk midplane region that contains the bulk of the disk 
mass and does not receive any direct UV illumination from the 
central star. This minimises further the effects of photodissocia- 
tion on our results. Therefore, we find it appropriate to use the 
12 CO mass derived from our best fit to the C 18 and 13 CO data 
to make an estimate of the total gas mass. 
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Table 1. The table shows the column of CO traced by 90% of the emission of each line, and the corresponding average temperature, 
at two radial intervals beyond 100 AU. These values are representative of the regions shown in Fig. [6] The column densities of 12 CO 
contained above the vertical locations z from Fig. [6] are shown. The average temperature values shown are column- weighted (i.e. 
mass-weighted) values over all cells which are located above the reported height z. 



Molecular 


z(AU) 


z(AU) 


N co (cm 2 ) 


Fraction Of 


N co (cnr 2 ) 


Fraction Of 


T ave (K) 


T ave (K) 


Line 


100-149 AU 


149-235 AU 


100-149 AU 


Total N co 


149-235 AU 


Total N co 


100-149 AU 


149-235 AU 


12 CO J=2-l 


12.7 


39.0 


6.5x10" 


7% 


6.6xl0 lb 


1% 


33.3 


38.6 


13 CO J=2-l 


7.6 


13.0 


4.2xl0 18 


17% 


2.6xl0 18 


18% 


29.1 


25.6 


C ls O J=2-l 


-2.5 


-4.3 


1.5xl0 19 


62% 


8.9xl0 18 


62% 


25.7 


21.8 



The abundance of 12 CO is a key factor in this con- 
version and the values found in the liter ature for the 
inters tellar medium range from 9.5x l0~ 4 dFrerking et al.l 
1982) to 2.0xl0~ 4 dLacv et all 119941) . If we adopt the 
12 CO abundance of 10~ 4 with respect to H2, we obtain 
M gas = 1.2 X M H , = (2.0^'g) x 10~ 2 M o . The error estimate in- 
cludes the error due to the rms of our C ls O data (41%) and 
to the isotopic ratio of [ 16 0]/[ 18 0] (4%). The upper margin is 
extended to include the 8% uncertainty in the gas mass esti- 
mate due to freeze-out. The factor 1 .2 takes into account the gas 
mass contained in helium, the second most important gas mass 
constituent after molecular hydrogen. Adopting a 12 CO abun- 
dance of 2.0xl0~ 4 with respect to H2 results in M gas as low 
as 6xl0~ 3 M . Therefore, due to the uncertainty largely domi- 
nated by the 12 CO abundance, we can only conclude that the gas 
mass of the disk around HD169142 falls in the range of (0.6- 
3.0)xl0- 2 M o . 



4.4. Gas to dust ratio 



In order to investigate the gas to dust ratio, we ran several 
models in which we assume M dllst =2.16xl0- 4 Mo, [CO]=10- 4 , 
12 C/ 13 C=77, 16 O/ 18 O=560 and gas to dust mass ratios of 100, 
50, 33 and 25. From the comparison of the synthetic visibili- 
ties of these models with the observations (Fig. [7] middle and 
lower panel) we find that a good fit for 13 CO and C ls O data is 
obtained for f g /<j=100. Under the assumptions used above, lower 
values of f g /d fail to provide a fit to the data. The 12 CO data 
(Fig. El top panel) shows a discrepancy, already discussed in sec- 
tion 4.3.1, with respect to the 13 CO and C 18 data in the sense 
that the model we use overestimates the 12 CO emission. This 
suggests that the model assumptions like the temperature, micro 
turbulence and 12 CO abundance used for the upper disk layers 
where the 12 CO emission arises from should be reconsidered. 
From the gas and dust masses derived in sections 4.2 and 4.3 we 
see that a gas to dust mass ratio of 100:1 only holds for special 
circumstances considered here: a maximum dust emissivity of 
2 cm 2 g du 1 st at 1.3 mm and a 12 CO abundance of 10~ 4 with re- 
spect to H2, or lower. For equally likely values of 1 cm 2 glL for 
the dust emissivity and 2xl0~ 4 for the 12 CO abundance, a gas to 
dust mass ratio of 25 would be inferred, suggesting that the disk 
would have lost significant amounts of gas with respect to dust. 
In fact, gas loss is expected in circumstellar disks due to photo- 
evaporation of gas by the stel lar radiation. The typical timescale 
for th is process of a few Myr ( Hol lenbach et al. 1994; Font et alj 
120041) is similar to the estimated age of HD169142 of 6+3 Myr 
dGradv et alJl2007h . 
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Fig. 7. Comparison of the observed vector-averaged 12 CO (top 
panel), 13 CO [middle panel) and C I8 (lower panel) line fluxes 
(black) to four models with different gas to dust ratios (red lines). 
The error bars represent the variance within each annular aver- 
age. The models assume the dust mass of 2.16xl0~ 4 M and gas 
to dust mass ratio f g /d=100,50,33,25 marked respectively by the 
full, dashed, dashed-dotted and dotted red lines in each panel. 
The line fluxes are integrated over the width of sj2.7 km s _1 (10 
spectral channels) centered on each line. 



4.5. Micro turbulence 

As mentioned in the previous section, the overestimate of the 
12 CO emission by the model which fits well the 13 CO and C ls O 



O. Panic et al.: Gas and dust mass in the disk around the Herbig Ae star HD169142 



9 



7= 2-1 line emission, may be due to model assumptions used 
for the upper disk layers, e.g., temperature, 12 CO relative abun- 
dance, turbulence. 

In section 4.1 we mention that a disk model with a less 
luminous A6 sp ectral type star also reproduces the SED of 
HD 169 142 well dDent et al.l l2006!). However, the temperature in 
that model is less than 10% lower than the temperature used in 
our model. This is insufficient to provide a significant decrease 
in the 12 CO line flux. 

A fit to all observed lines simultaneously by varying the iso- 
topic ratios, yields an unusualy low value of both the 12 C/ 13 C 
and the 16 0/ 18 isotopic ratio: 21 and 171, respectivelly. These 
values are far from the observationally derived ratios in the in- 
terstellar medium and this scenario would be highly unlikely. 
An over-abundance of 13 CO was measured at large radial dis- 
tances from the star (several hundreds of AU) in the disks DM 
Tau, MWC 480 and LkCa 15 dPietu et alj2007l) In general, disks 
are often found to have 13 CO and C 18 abundance s depleted by 
large factors , mainly in disks around T Tauri stars dDutrev et alj 
I1994T 119961) . 

In the above analysis, the turbulent line FWHM was set to 
0. 16 km s _I throughout the disk. In a more realistic scenario, the 
turbulence in the cold and dense midplane differs from the tur- 
bulence in the warmer and less dense surface layers. We investi- 
gate the effect of micro turbulence on the modeled line emission 
and find that a simultaneous match is possible in our model, as 
shown in Fig. [8] This fit is obtained for a gas to dust ratio of 
33, i.e., a disk gas mass of 7xl0~ 3 M , if a higher value of tur- 
bulence is assumed for the disk midplane (turbulent FWHM of 
0.33 km s~') and no turbulence for the upper disk layers traced 
by the 12 CO emission (0.00 km s _1 , only thermal broadening). 
For the 12 CO J=2-l line, that is fully optically thick at all radii 
in our model, the turbulent line broadening is degenerate with 
temperature and density in its effect on the line flux and profile. 
This is because the turbulence additionally broadens the ther- 
mally broadened line profile and in this way the optical thick- 
ness decreases allowing for more gas column to be traced. On 
the other hand, by decreasing the turbulence, at a given tempera- 
ture, the same 12 CO data can be fit with comparably more Mco- 
The effect of the increase in turbulence on the 13 CO and C ls O 
lines is present but less pronounced, due to their lower optical 
thickness. In particular, the C 18 7=2-1 line emission is only 
marginally optically thick. The effect of turbulence on the mod- 
eled C ls O .7=2-1 line emission is below the rms of our obser- 
vations. It is not possible to obtain simultaneous fits to all data 
if the gas to dust ratio is higher than 33 (disk mass higher than 
7xl0~ 3 M ). A lower gas to dust ratio (disk mass lower than 
7xl0~ 3 M ) would require the turbulent FWHM of the lines to 
be larger than 0.33 km s _1 in disk midplane. 

5. Summary 

The circumstellar disk around the Herbig Ae star HD 1 69 142 was 
imaged in 13 CO and C 18 7=2-1 line emission at 1'.'4 resolution. 
The disk is resolved and the emission extends to 200-300 AU 
from the star at 2<x levels. We adopt an accretion disk model 
which matches the SED of HD169142 and we use a molecular 
excitation and radiative transfer code to analyse our data. We re- 
visit existing 1.3 mm and 12 CO data and interpret it in the view 
of our new results. The observed kinematic pattern and line pro- 
files are consistent with a disk seen at 13° inclination from face 
on, in Keplerian rotation around a 2 M star. The structure of the 
emission suggests a disk radius of roughly 200 AU for both gas 
and dust. 




60 80 too 
UV distance (kA) 

Fig. 8. Comparison of the observed vector-averag ed 12 CO (top 
panel), 13 CO (middle panel) and C 18 (lower panel) line fluxes 
(black) to a model which assumes the total disk mass of 
7xl0~ 3 M , gas to dust mass ratio f g /d=33, no microturbulent 
line broadening in the upper disk layers from which the 12 CO 
emission originates (only thermal broadening), and a microtur- 
bulent broadening of 0.33 km s -1 FWHM in the layers from 
which the I3 CO and C 18 emission arises. 



We calculate the 



of 



minimum disk dust mass 
M dust =2.16xl0~ 4 M through the fit to the 1.3 mm 
data using the maximum dust opacity theoretically predicted 
in circumstellar disks. 

The adopted disk model reproduces our molecu- 
lar line observations with M c is o =4.6xl0~ 8 M and 
Mi3 CO =2.9xlO" 7 M . The derived amount of C 18 implies 
a M co =(2.4+1.0)xl0- 5 M . The 13 CO amount is also in 
agreement with this value. 

We identify the regions of the disk that contribute most to the 
emission of I2 CO, 13 CO and C ls O J=2-l line and find, as ex- 
pected due to different optical depth of these lines, that they 
are physically distinct. In particular, the C 18 J=2-l line is 
marginally optically thin throughout the disk and traces the 
bulk of the mass contained in the disk midplane. 
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- We find that the effect of freeze-out in the model is relatively 
small, allowing us to make a gas mass estimate based on our 
observations. For plausible CO abundances of (l-2)xl0~ 4 , 
we can only derive a range of M gas =(0.6-3.0)xl0~ 2 M Q . 

- The comparison of the derived gas and dust mass shows that 
the gas to dust mass ratio of 100 is only possible under as- 
sumption of a high dust opacity of 2 cm 2 g di J jr and low 12 CO 
abundance of 10" 4 with respect to H2. Otherwise, this ratio 
may be as low as 25. 

- Existing observations of the optically thick 12 CO line can be 
fit simultaneously to our 13 CO and C ls O data by the model 
adopted here, only if the micro turbulence in the disk mid- 
plane is higher than the micro turbulence in the upper disk 
layers which dominate the 12 CO emission. 

The discrepancy between the CO gas mass derived here and 
that derived based on the 12 CO line emission likely implies that 
the micro turbulence is much lower than 0.1 km s _1 in the up- 
per disk layers where the 12 CO emission originates. Better con- 
straints of the 12 CO abundance in disks are needed in order to 
enable us to measure gas mass more precisely in disks. Better 
knowledge of dust properties such as dust size distribution and 
composition is essential for a better estimate of dust opacity in 
disks. However, these properties vary from source to source and 
are not necessarily linked with other parameters, age or mass 
for example and the studies need to focus on individual sources. 
Since the disks around Herbig Ae stars have the advantage that 
the gas mass estimate from C 18 line emission is more reliable, 
it would be useful to investigate dust properties, i.e., emissivity 
in these sources in more detail to be able to have a reliable esti- 
mate of dust mass as well. This would, for the first time, open a 
real opportunity to study the gas to dust mass ratio, the unknown 
property on which the theories of disk evolution and planet for- 
mation depend strongly. 
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